Entanglement in spin-one Heisenberg chains 
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By using the concept of negativity, we study entanglement in spin-one Heisenberg chains. Both 
the bilinear chain and the bilinear-biquadratic chain are considered. Due to the SU(2) symmetry, the 
negativity can be determined by two correlators, which greatly facilitate the study of entanglement 
properties. Analytical results of negativity are obtained in the bilinear model up to four spins 
and the two-spin bilinear-biquadratic model, and numerical results of negativity are presented. We 
determine the threshold temperature before which the thermal state is doomed to be entangled. 
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I. INTRODUCTION 

Since Haldanc predicts that the one-dimensional 
Heisenberg chain has a spin gap for integer spins 0, 
the physics of quantum spin chains has been the subject 
of many theoretical and experimental studies. In these 
studies, the bilinear spin-one Heisenberg model and the 
bilinear-biquadratic Heisenberg model have played im- 
portant roles |2j, y, |4| . The corresponding Hamiltonians 
are given by 

JV 



H 1 =Y,JSi-S i+ U (!) 

N 

#2 =J2 [ JS > ' S '+ 1 + ^ ■ Sl + 1 ) 2 ] » ( 2 ) 



respectively. Here, we have assumed the periodic bound- 
ary condition, and obviously, these two Hamiltonians 
exhibits a SU(2) symmetry. Moreover, the bilinear- 
biquadratic model exhibits very rich phase diagram 

Recently, the study of entanglement properties in 
Heisenberg systems have received much attention |r3|-[39|. 
Quantum entanglement lies at the heart of quantum me- 
chanics, and can be exploited to accomplish some phys- 
ical tasks such as quantum teleportation |40| . Spin-half 
systems have been considered in most of these studies. 
However, due to the lack of entanglement measure for 
higher spin systems, the entanglement in higher spin sys- 
tems have been less studied. There are several proceed- 
ing works on entanglement i n sp in-one chains. Fan et 
al. [33| and Verstraete et al. |3J] studied entanglement 
in the bilinear-biquadratic model with a special value of 
7 = 1/3, i.e., the AKLT model @. Zhou et al. studied 
entanglement in the Hamiltonian H2 for the case of two 
spins |35| . 

In this paper, by using the concept of negativity |4l| . 
we study pairwise entanglement in both the bilinear and 
the bilinear-biquadratic Heisenberg spin-one models. For 
the case of higher spins, a non-entangled state has nec- 
essarily a positive partial transpose (PPT) according to 
the Peres- Horodecki criterion 1431 . In the case of two 



spin halves, and the case of (1/2,1) mixed spins, a PPT 
is also sufficient. However, in the case of two spin-one 
particles, a PPT is not sufficient. Nevertheless, the neg- 
ative partial transpose (NPT) gives a sufficient condition 
for entanglement, and due to the SU(2) symmetry in the 
systems, the NPT is expected to fully capture the entan- 
glement properties. 

The Peres-Horodecki criterion give a qualitative way 
for judging if the state is entangled. The quantita- 
tive version of the criterion was developed by Vidal and 
Werner |4l|. They presented a measure of entanglement 
called negativity that can be computed efficiently, and 
the negativity does not increase under local manupula- 
tions of the system. The negativity of a state p is defined 
as 



Mi 



(3) 



where pi is the negative eigenvalue of p T2 , and T2 denotes 
the partial transpose with respect to the second system. 
The negativity Af is related to the trace norm of p Ta via 



Af{ P ) 



\\p T *h 



(4) 



where the trace norm of p 1 " 2 is equal to the sum of the 
absolute values of the eigenvalues of p T2 . If j\f > 0, then 
the two-spin state is entangled. 

We study entanglement in both the ground state and 
the thermal state. The state of a system at thermal 
equilibrium described by the density operator p{T) — 
exp(—f3H)/Z, where f3 — 1/fegT, fee is the Boltzmann's 
constant, which is assume to be 1 throughout the pa- 
per , and Z = Tr{exp(— (3H)} is the partition function. 
The entanglement in the thermal state is called thermal 
entanglement. 

We organize the paper as follows. In Sec. II, we give 
the exact forms of the negativity for an SU(2)-invariant 
state, and show how the negativity is related to two cor- 
relators. We also give that how to obtain negativity from 
the ground-state energy and partition function in the 
bilinear-biquadratic model. We study entanglement in 
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the bilinear and bilinear-biquadratic models in Sec. Ill 
and IV, respectively. Some analytical and numerical re- 
sults of negativity are obtained. We conclude in Sec. V. 



II. NEGATIVITY AND CORRELATORS 

Schliemann considered the entanglement of two spin- 
one particles via the Peres- Horodecki criteria [3l|. and 
find that the SU(2)-invariant two-spin state is entangled 
if either of the following inequalities holds 



((Si • S,) 2 ) > 2, 
((S 4 -S,) 2 ) + (S 4 -S,) < i. 



(5) 



Now we explicitly give the expression of negativity for 
the SU(2)-invariant two-spin state. 

According to the SU(2)-invariant symmetry, any state 
of two spin-one particles have the general form 

p =G\S = Q,S Z = 0)(S = 0,S Z = 0| 



H 



]T \s = i,s z ){s = i,s z 



s.=-i 



1-G-H 



\S = 2,S Z ){S = 2,S Z \, (6) 



Sz=-2 



where \S, S z ) denotes a state of total spin S and z com- 
ponent S z , and 



G=i[((S,.S J ) 2 )-l], 



(7) 



In order to perform partial transpose, the product ba- 
sis spanned by {ISi = 1, S\ z } (g> \S2 = 1, S2 Z )} is a natural 
choice. By using the Clebsch-Gordan coefficients, we may 
write state p in the product basis. The partially trans- 
posed with respect to the second spin p T2 can be written 
in a block-diagonal form with two lxl block, two 2x2 
block, and one 3x3 block. After diagonalization of each 
block, one find that the following only two eigenvalues of 
p T2 are possibly negative [3l|. 



Mi=i(2-((S l .S,) 2 )), 
Ai2 =i((S 4 .S,) + ((S l .S J ) 2 )-l) 



(8) 



Moreover, p\ and fi2 occur with multiplicities 3 and 1, 
respectively. Therefore, the negativity is obtained as 

A^> =imax[0,<(S J -S J ) 2 )-2] 

+i max[0, 1 - ((Sj • Sj)) - ((Si • S,) 2 )]. (9) 



We see that for the SU(2)-invariant state, the negativity 
is completely determined by two correlators ((Sj ■ Sj)} 
and ((S.-S,) 2 ). 

Recall that the swap operator between two spin-one 
particles is given by 



&ij — Si • Sj 



(Si 



s,) 2 



(10) 



where / denotes the 9x9 identity matrix. Then, the 
negativity can be written in the following form 



AT™ =imax[0,(5 y )-(S. i -S,)-l] 
+- max[0, -(«%)]. 



(11) 



We see that if the expectation value (Sij) < 0, the state is 
entangled. The swap operator satisfies Sy = 1, and thus 
it has only two eigenvalues ±1. If a state is a eigenstate of 
the swap operator, the expression l|ll|) can be simplified. 
When the corresponding eigenvalue is 1, Equation l|ll|) 
simplifies to 



= imax[0,-(Si • S 



(12) 



and when the eigenvalue is -1, the equation simplifies to 
1 1 



3 



max 0, 



(Si ' 



2]- 



(13) 



In the former case, the state is entangled if (Sj • Sj) < 0, 
and in the latter case, the state is an entangled state, and 
the negativity is larger than or equal to 1/3. 

Now we consider the bilinear-biquadratic spin-one 
Heisenberg model described by the Hamiltonian Hi- By 
applying the Hellmann-Feynman theorem to the ground 
state of H2 and considering the translational invariance, 
we may obtain the correlators as 



GS 



((8,-Si+i)) - ^^fp, ((S s ;-S l+1 ) 2 > - 



(14) 



where Eqs is the ground-state energy. Substituting the 
above equation into Eq. © yields 



0. 



1 dE t 



GS 



0,1 



N d 7 

1 dE GS 



1 dE, 



GS 



N dJ 



N 9 7 



(15) 



For the case of finite temperature, we have 



0. 



-1 dZ 



N(3Z <9 7 

1 dz 

0,1 + 



1 dZ 



NfiZ dJ N(3Z 97 



(16) 



We see that the knowledge of ground-state energy (par- 
tition function) is sufficient to determine the negativity 
for the case of zero temperature (finite temperature). 
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III. BILINEAR HEISENBERG MODEL 

Let us know consider the entanglement in the bilinear 
Heisenberg model. Due to the nearest-neighbor character 
of the interaction, the entanglement between two nearest- 
neighbor spins is prominent compared with two non- 
nearest-neighbor spins. Thus, we focus on the nearest- 
neighbor case in the following discussions of entangle- 
ment. 



A. Two spins 

For systems with a few spin, we aim at obtaining an- 
alytical results of negativity. The Hamiltonian for two 
spins can be written as 

Hi = S 1 -S 2 -i[(S 1 + S 2 ) 2 -S2-S2], (17) 

from which all the eigenvalues of the system are given by 



E = -2(1), Ej. = -1(3), E 2 = 1(5), 



(18) 



where the number in the bracket denotes the degeneracy. 

We investigate the entanglement of all eigenstates of 
the system. When an energy level of our system is non- 
degenerate, the corresponding eigenstate is pure. When 
a fc-th energy level is degenerate, we assume that the 
corresponding state is an equal mixture of all eigenstates 
with energy Thus, the state correspoding to the fc- 
th level with degeneracy becomes a mixed other than 
pure, keeping all symmetries of the Hamiltonian. A 
degenerate ground state is called thermal ground state 
in the sense that it can be obtained from the thermal 
state cxp[— H / (ksT)]/ Z by taking the zero-temperature 
limit [15|. The fc-th eigenstate can be considered as 
the thermal ground state of the nonlinear Hamiltonian 
H' given by W = (H - E k ) 2 . Note that Hamiltonian H' 
inherits all symmetries of Hamiltonian H. 

As we consider interaction of two spins, from Eqs. lfT7l) 
and we obtain another form of the negativity as 

N {12) =imax[0,(i/ 1 2 ) -2] 

+ imax[0,l-Cff 1 +tf 1 2 )]. (19) 



Having known negativities of all eigenstates, we next 
consider the case of finite temperature. The cumu- 
lants can be obtained from the partition function. From 
Eq. i|18|) . the partition function is given by 



Z = e 2 ? 



(21) 



A cumulant of arbitrary order can be calculated from the 
partition function, 



(-1)™ d n 



Z d(3 n 
(-1) 



Z 



z 



[2"e 2/3 + 3e' 3 + 5(-l)"e-' 3 ] 



(22) 



Substituting the cumulants with n = 1,2 to Eq. Ijl9|l 
yields 



TV max (0, 2e 20 - 3e^ - be~ f) ) 



— max (0, 3e p - e 2/3 



5e" 



(23) 



Thus, we obtain the analytical expression of the negativ- 
ity. 

The second term in Eq. (|23|) can be shown to be zero. 
To see this fact, it is sufficient to show that F(x) — x 3 — 
3x 2 + 5 > 0, where x = e 13 > 1. It is direct to check 
that the function F takes its minimum 1 at x = 2. As 
the minimum is large than zero, the function is positive 
definite. Thus, equation (|23|l simplifies to 



1 

2Z 



max (0, 2e 2/3 - 3e p - 5e 



(24) 



The behavior of the negativity versus temperature is sim- 
ilar to that of the concurrence ^3 m the spin- half Heisen- 
berg model Q, namely, the negativity decreases as the 
temperature increases, and there exists a threshold value 
of temperature T t h, after which the negativity vanished. 
This behavior is easy to understand as the increase of 
temperature leads to the increase of probability of the 
excited states in the the thermal state, and the excited 
states are less entangled in comparison with the ground 
state. From Eq. Ij24(l . the threshold temperature can be 
analytically obtained as 



T th =- 



1 



ln(i 



-2 1 2(ll+2V30) 1 /3 

d. 3667. 



(ll+2y / 30) 1 /3 n 
2 > 



(25) 



To determine the negativity, it is sufficient to know the 
cumulants (Hi) and {H 2 ). 

From Eqs. I|19l) and l|18fl . the negativities correspond 
to the fc-th level are obtained as 

A^ 12) = 1, Af[ 12) - 1/3, A/f 2) - 0. (20) 

We see that the ground state is a maximally entangled 
state, the first-excited state is also entangled, but the 
negativity of the second-excited state is zero. 



B. Three spins 



The Hamiltonian for three spins is rewritten as 

H l = i[(Si + S 2 + S 3 ) 2 - S 2 - S 2 - S§], (26) 

from which the ground-state energy and the correlator 
(Si • S2) are immediately obtained as 



GS 



-3, (Si ■ S 2 



1. 



(27) 
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In order to know the ground-state negativity, we need to 
calculator another correlator ((Si • S 2 ) 2 ). 

By considering the translational invariance and using 
similar techniques given by Refs. [tH liH |4^| . the ground- 
state vector is obtained as 



1 



Wgs = ^(|012) + |201) + |120) 
-|021> - |102> - |210», 



(28) 



where \n) denote the state \s = l,m = s — n), the com- 
mon eigenstate of S 2 and S z . Then, we can check that 

Si • S 2 |*) GS = -|*)gs- (29) 

Thus, the correlator ((Si • S2) 2 ) is found to be 

((Si • S 2 ) 2 ) = 1. (30) 

Substituting Eqs. J23 and (EOJ to Eq. @ yields 

AT =1/3. (31) 

We see that spins 1 and 2 are in an entangled state at 
zero temperature. With the increase of temperature, the 
negativity monotonically decreases until it reaches the 
threshold value T t h — 0.9085, after which the negativity 
vanishes. 

C. Four spins 

Now we consider the four-spin case, and the corre- 
sponding Hamiltonian can be written as 

Hx = l[(Si + S 2 + S 3 + S 4 ) 2 - (Si + S3) 2 - (S 2 + S 4 ) 2 ]. 

The standard angular momentum coupling theory di- 
rectly yields the ground-state energy and the correlator 
(Si • S 2 ) 



E, 



GS 



-6, (Si -S 2 



-3/2. 



(33) 



Then, we need to compute another correlator ((Si • S 2 ) 2 ) 
or alternatively the expectation value (£12)- So, it is 
necessary to know the exact form of the ground state. 

By using similar techniques given by Refs. |4Hl45ll4^ . 
the ground-state vector is obtained as 



|¥>gs =1/2|Vi> - 3/2|^ 2 > + |V>3) 
- 3/2^4} + 3/V2|V> 5 ) + 



|V>6 



(34) 



where 





=1/2(|0022) + 


|2002) - 


h |2200) - 


- |0220)), 


1^2) 


=1/2(|0112) + 


|2011) - 


h |1201) - 


-11120)), 


1^3) 


=1/2(|0121) + 


|1012) - 


h |2101) - 


-11210)), 




=1/2(|0211) + 


|1021) - 


h |1102) - 


-12110)), 




=1/V2(|0202) 


+ |2020)), 






= |1111). 










FIG. 1: Negativity versus 7 for different temperatures in the 
two-spin model ( J = 1). 



Then, from the explicit form of the ground state, after 
two-page calculations, we obtain the expectation value of 
the swap operator as 

(S12) = 1/6. (36) 

Substituting Eqs. and (JSgJl to Eq. (JTTJ) leads to 

A/ - =1/3. (37) 

It is interesting to see that the ground-state negativity 
in the four-qubit model is the same as that in the three- 
qubit model. The threshold value can be found to be 
Tth = 1.3804. 

For N > 5, it is hard to obtain analytical results of neg- 
ativity. The behaviors of negativity are similar to those 
for N < 4, namely, with the increase of temperature, the 
negativity decreases until it vanishes at threshold tem- 
perature Tth- For instance, the threshold temperatures 
T t h ~ 0.95 and T t h ~ 1.21 for five and six spins, respec- 
tively. The negativity for two nearest-neighbors spins is 
estimated as TV = 0.1240 (W = 0.2509) for the case of 
five spins (six spins). 



IV. BILINEAR-BIQUADRATIC SPIN-ONE 
HEISENBERG CHAIN 



We now study entanglement properties in the bilinear- 
biquadratic spin-one Heisenberg model, and first consider 
the case of two spins. From Eq. (|15|) with N — 1, if we 
know the ground-state energy, the negativity is readily 
obtained. The ground-state energy is given by 



Egs = 



2 J 
-1J 



■ 47 
7 



when 7 < 1/3, 
when 7 > 1/3, 



(38) 



(35) 



We see that there exits a level crossing at the point of 
7 = 1/3. Then, substituting the above equation into 



FIG. 2: Negativity versus 7 for different temperatures in the 
three-spin model (J = 1). 



FIG. 3: Negativity versus 7 for different temperatures in the 
four-spin model ( J = 1). 



Eq. iJTSJ yields 



TV 



1 when 7 < 1/3, 
1/3 when 7 > 1/3, 



(39) 



Before the point 7 = 1/3, the negativity of the ground- 
state is 1, while the negativity of the first-excited state is 
1/3. After the cross point, the ground and first-excited 
interchanges, and thus, the negativity of the ground state 
after the cross point is 1/3. It is interesting to see that 
the model at the cross point is just the AKLT model. 

In Fig. 1, we plot the negativity versus 7 for different 
temperatures. The level cross greatly affects the behav- 
iors of the negativity at finite temperatures. For a small 
temperature (T = 0.05), the negativity displays a jump 
from 1 to 1/3 near the cross point. For higher temper- 
atures, the negativity first decreases, and then increases 
at 7 increases from -1 to 1. For T = 1.5, we observe that 
there exists a range of 7, in which the negativity is zero. 

For the three-spin case, we plot the negativity versus 
7 for different temperatures in Fig. 2. For a low temper- 
ature T = 0.015, we observe a dip, which results from 
the level crossing near the point of 7 = —0.2121. When 
T = 0.1, the dip becomes more evident. For the cases of 
higher temperatures (T — 0.5 and T — 1.0), there exists 
a range of parameter 7, in which the negativity is zero. 

For the four-spin case, we also a plot of the negativity 
for different temperatures. For T = 0.03, as 7 increases, 
the negativity decreases until it reaches its minimum, and 
then increases. For T = 0.5 and 1.0, the behaviors of 
negativity are similar to the case of T = 0.01, and the 
difference is that the minima shifts left. There are some 
common features in the behaviors of negativity for differ- 
ent number of spins. The maximum value of negativity 
occurs at 7 = —1; for higher temperatures, there exists 
a range of 7, in which the negativity is zero. 

From Figs. 1-3, we observe that the thermal state is 
always entangled at a lower temperature. When temper- 



ature increases, the negativity decreases until it reaches 
zero, namely, the thermal fluctuation suppresses entan- 
glement. Before the threshold temperature, the state 
is doomed to be entangled. We numerically calculated 
the threshold temperature and the result are shown in 
Fig. 4. The threshold temperature decreases nearly lin- 
early when 7 increases from -1 to a certain value of 7. 
After reaching a minimum, it begin to increase. We see 
that the behaviors of the threshold temperature are sim- 
ilar for different number of spins. 

As a final remark, we consider the following Hamilto- 
nian 



n 1 n y 

1+3 \i=\ / 



(40) 



where the interaction is between all spins, and there are 
all together N(N — l)/2 terms. The system not only 
shows a SU(2) symmetry, but also an exchange sym- 
metry, namely, the Hamiltonian in invariant under ex- 
change operation SijH 3 Sij — H 3 . For N — 2,3, the 
model is identical to Hamiltonian H2 ■ We know that the 
ground state is non-degenerate when N = 2,3, and thus 
it must be an eigenstate of £y and Eqs. (|12fl and i|13|) 
can apply. From the angular momentum coupling the- 
ory, the ground-state energy of H3 is readily obtained as 
E GS = -N, and thus we have (S, ■ Sj) = -2/(N - 1). 
Then, from Eqs. (|12|) and l|13fl . the negativity can be ei- 
ther 1/(N - 1) or 1/3. For N = 2 (N = 3), the ground 
state is symmetric (antisymmetric) and then the negativ- 
ity is 1 (1/3), consistent with previous results. However, 
for N > 4, the ground-state is degenerate and we cannot 
apply Eqs. (|12|) and (|13fl . The numerical results show 
that the negativity is zero for N > 4. 
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y 



FIG. 4: Threshold temperature versus 7 for different number 
of spins. 

V. CONCLUSIONS 

In conclusion, by using the concept of negativity, we 
have studied entanglement in spin-one Heisenberg chains. 
Both the bilinear model and bilinear-biquadratic model 
are considered. Although NPT only give a sufficient con- 
dition for entanglement, due to the SU(2) symmetry, we 
believe that this condition considerably captures entan- 
glement properties of the system. Moreover, the nega- 



tivity gives an upper bound to the distillation of entan- 
glement 1 1 lj , one of the fundamental entanglement mea- 
sures. We have given explicitly the relation between the 
negativity and two correlators. The merit of this relation 
is that the two correlators completely determine the neg- 
ativity and it facilitates our discussions of entanglement 
properties. 

We have obtained analytical results of negativity in 
the bilinear model up to four spins and in the two- 
spin bilinear-biquadratic model. We numerically calcu- 
lated entanglement in the bilinear-biquadratic model for 
N = 2,3,4, and the threshold temperatures versus 7 
are also given. We have restricted us to the small-size 
systems, and aimed at obtaining analytical results via 
symmetry considerations and getting some numerical re- 
sults via the exact diagonalization method. However, 
for larger systems, the exact diagonalization method is 
not a viable route. It is interesting to investigate large 
systems by some mature numerical methods such as the 
quantum monte-carlo method and density-matrix reno- 
malization group method. And it is also interesting to 
consider other SU(2)-invariant spin-one systems such as 
the dimerized and frustrated systems. 
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